%plot Inertia of gait spline
tmax = time(6);
the2 = PP;
for t=0:0.001:tmax
    [Th, Tk, Jh(round(t/0.001)+1), Jk(round(t/0.001)+1)] = load_analysis(0, the2(round(t/0.001)+1));
end
t = 0:0.001:tmax;
plot(t, Jh.*1e7)